Space-time properties of free motion time-of-arrival eigenstates. 
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The properties of the time-of-arrival operator for free motion introduced by Aharonov and Bohm 
and of its self-adjoint variants are studied. The domains of applicability of the different approaches 
are clarified. It is shown that the arrival time of the eigenstates is not sharply defined. However, 
strongly peaked real-space (normalized) wave packets constructed with narrow Gaussian envelopes 
centred on one of the eigenstates provide an arbitrarily sharp arrival time. 
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I. INTRODUCTION 



The extension of the classical arrival time concept to quantum mechanics is problematic because of the absence 
of trajectories in the standard interpretation of the quantum formalism. The question is important since detection 
of particles in time-of-flight and coincidence experiments are quite common, and quantum mechanics should be able 
to predict and interpret the statistics of the arrivals. A recent review discusses the main theoretical approaches and 
open questions jlj]. 

The customary way to specify a dynamical variable in quantum mechanics is to define the corresponding operator, 
but Pauli showed that a self adjoint time operator conjugate to the Hamiltonian H cannot be constructed if the 
spectrum of H is bounded from below §, see also There are however several ways to circumvent this objection 

and quantize the classical arrival time, at least in simple cases such as free motion. Aharonov and Bohm introduced 
the operator || 
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T = MH--^)=--\^ + -ji) , (i) 



by symmetrizing the classical expression for the time-of-arrival at X — for a particle with position q and momentum 
p at time t = 0. The same result is obtained by using the quantization rules of Weyl, Rivier or Born- Jordan [Q, 
but this coincidence does not justify ([!]) since, a priori, no rule is more fundamental than any other rule. They 
simply provide a number of quantum generalizations of a classical quantity; one (or some) of these can be selected by 
requiring either agreement with a particular experimental procedure or that certain desirable conditions be satisfied, 
such as properties obeyed by the classical counterpart. It is thus essential to carefully examine the behaviour of a 
given operator and of its eigenstates to ascertain their ultimate physical content. Paul and others have examined 
the mathematical properties of T as a linear operator in Hilbert space In section II we shall first review 

these properties. Then the relations between T and positive operator valued measures, or its self-adjoint variants are 
examined, and the important question of the domains of applicability of the theories associated with the different 
operators is clarified. In section III less abstract aspects are studied, such as the functional form of the eigenstates in 
the coordinate representation and their time dependence. Finally, the behaviour of normalized wave packets "peaked" 
around one of the eigenstates is discussed in section IV. 



II. TIME OPERATORS 

T is not a self-adjoint operator but it is in a sense the next best thing, namely, a maximal symmetric operator 
@ fj||. This means that it is hermitean (or "symmetric" in the mathematical literature) but it does not admit a 
self-adjoint extension. Indeed a first way to circumvent Pauli's argument is to relax the standard requirement and 
admit hermitean but non self-adjoint operators as representations of physical properties. By imposing hermiticity 
and square integrability of the states in its range, its domain D(T) is given by the set of states with momentum 
representation obeying || 

(2) 
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The eigenvalue equation in momentum representation, {p\T\T) = T(jp\T), 

'I d 1 
^p dp 2p 2 

is solved (for p ^ 0) by 

(p\T, a) = (|p|/m^) 1 / 2 e lp2T /( 2m?l )e(ap) , (4) 

with a = ±. These (improper) eigenfunctions have the appropriate transformation behaviour for a state arriving at 
time T, 

e -&*/*\T,a) = \T-t,a), (5) 

namely, if the original state had already evolved for a time t then it would arrive at T — t instead; but they are not 
square integrable functions so they do not define physically realizable state vectors in the Hilbert space Ti, = L 2 . The 
constant factor in (Q) is chosen to satisfy the resolution of the identity 

/oo 
dT\T,a)(T,a\. (6) 

In the physical interpretation of this expression, one assumes that all particles arrive sooner or later. Notice the 
inclusion of negative times. In an experimental context, all measured arrival times are positive if the single particle 
state ip(t) of interest is prepared at the instant t = 0. However, as discussed in considerable detail by Grot et al. Jl3| , 
the theoretical analysis is much simpler if one considers a different problem in which it is imagined that the state is 
prepared at t — — oo in the state which (in the absence of interaction) evolves to the actually prepared state ip(t = 0), 
and the arrival time distribution includes both negative as well as positive arrival times. This is very useful from the 
mathematical point of view (e.g. see Eq. (^), or Eqs. ( fi"5| ) and ( p2| ) below.) 
The eigenfunctions \T,a) are complete, but they are not orthogonal, 



\T , a. \T,a) = — — 



8{T - T') +-V ' 
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(7) 



The non orthogonality of the eigenstates has been associated with the intrinsically "unsharp" character of this "observ- 
able" by the proponents of the "operational approach to quantum mechanics" [0 . These authors, and Gianitrappani 
in particular M, point out that the states \T, a) provide a "positive operator valued measure" (POVM) for the arrival 
time. This means that for an interval of time [Ti,Ta] the positive, bounded operator 

B(T 2 ,T 1 ) = yf 2 \T,a)dT{T,a\ (8) 



Ti 



can be constructed. This is not a proj ector, but one does not necessarily require projectors in order to introduce 
probabilities in quantum mechanics IHQi^I . By taking the trace with the normalized physical density operator p, the 
function 

P^T^^triB^T^p] (9) 

fulfills in principle the conditions of a probability for arrival between T± and it is positive, additive for disjoint 
sets, and, using (^|), tends to 1 as Ti — > — oo and T 2 — > oo. In particular, an arrival time distribution (for a pure state 
-0 at time t) is defined by 

Il(T^(t))=J2\(T,<x\m}\ 2 - (10) 

ot 

According to (j|), the POVM satisfies the "covariance condition" 

e -iHt/H^ T ^ Tx yHt/h = £( Ta _ t Tl _ t y (11) 

and this means in particular that 
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n(T;^(Q))=n(T-t;^(t)), 



(12) 



In words, the probability for arriving at T is equal to the probability for arriving at T — t when the original state has 
evolved a time t. This is a basic physical requirement that any candidate for an arrival time distribution should obey. 
This is of course also true in classical mechanics. 

The time-of-arrival operator T is the "first moment" of the POV measure. For any 4> £ H and ip E D(T), 

ee / T{<j>\B{dTm = (4\T\il)) =J2 I dT M r > «> T < T ' «l^> ■ ( l3 ) 

There are in principle other POV measures compatible with (|ll] ) and (|l3|), but (|§|) has the following unique feature: 
Its second moment operator obeys, for states in D(T), 

Aft/*) ee / T 2 (V,|%T)|^) - (TtP\TiP) = , (14) 



which is by no means an obvious relation since the eigenvectors of T are not orthogonal (An equation like (|lj) would 
be satisfied trivially by a self-adjoint operator since in that case the moments of the distribution of the observable 
are obtained as the expectation values of the powers of the operator.) A(-0) is called the "variance form" JhJ. In 
classical mechanics the true arrival time distribution makes the analogous quantity minimal so this property has also 
been invoked to select a proper quantum distribution 15 ic|]. 



Another condition satisfied by II (T) is that, for the state ipi defined by (p\ipi) = (p\ip)* , 

n(T,v(o)) = n(-T,^(o)), (15) 

which follows immediately from the symmetry property 

{p\T ) a)* = {p\-T,a). (16) 

Eq. (|l5]) is also a classically motivated relation that must hold when the arrival point is X — |L5|,|l0|. (If only 
covariance and minimum variance were imposed the arrival point would not be specified.) 

The distribution of arrival times (jlpj) , was first obtained by Kijowski for states with positive momenta or negative 
momenta o nly (a = 1 or a — — 1 but not both) |L5JL and has been later rederived, studied or generalized by several 



authors @pUT|J||3@- Werner in particular £0 
uniqueness of ( 



with the same restriction imposed by Kijowski, justified the 
( lot) subject to the conditions (|l2]), (fL4), and (jl^). Grot, Rovelli and Tate |l^] introduced a regularized 
self-adjoint operator and considered the full expression (|Io| ) for possible application to more general states having both 
positive and negative momenta but vanishing in the proximity of p — 0. Dclgado and Muga with the restriction 
imposed by Kijowski and Werner (i.e., for states with purely positive or negative momenta) arrive at the distribution 
using a different self-adjoint operator. We shall next discuss these two proposals for self-adjoint arrival time operators. 



A. SELF-ADJOINT VARIANTS OF T 

Grot, Rovelli and Tate Jl3| trace the non-orthogonality of the eigenstates to the singularity at p — 0, and produce 
a "regularized" self adjoint time operator T e with eigenstates that differ from (^]) only in a small momentum region 
around p = 0, 

(p\T±) £ = Q(±p){hmf £ (p)r 1/2 exp f , (17) 

where e is a small positive number, 



mh J± £ f e (p')J ' 



■^-{l-^p, \p\<t' (18) 

and O is the Heaviside "step" function. The integral in (|lj]), that plays the role of the energy when divided by m, 
becomes 

P J*L_ f(p 2 -e 2 )/2, |p|> £ {m 
±efe(p') ~ \e 2 H\p\/e), \p\<e> 
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so that the regularization amounts to changing the energy spectrum in the region around p = by introducing 
negative energies. In this manner Pauli's objection is avoided. While Grot, Rovelli and Tate modify the eigenstates 
(and therefore the time operator) in the proximity of p — 0, Paul discussed a similar idea ||, namely, to modify 
the physical state vector wave functions infinitcsimally around this point so that T can be applied to them. He 
concluded however that such modifications were not physically meaningful since these infinitesimal changes may lead, 
for example, to arbitrary values of the expectation value of the square of T, (T 2 ). Note that the small p region is 
responsible for the long time asymptotic behaviour of the arrival time distribution (see the Appendix A) and therefore 
any infinitesimal change there affects drastically quantities such as (T™), even though it will affect only infinitesimally 
other expectation values, for example (q n ), or (p n ), with n = 1,2,3... p|. Paul's observation is correct, but we shall 
argue in the discussion that the consequences are not necessarily as negative as he thought. 

A second proposal that was already pointed out by Kijowski (section 8 of fl5|| ) and has been developed further by 
Delgado and Muga || is to take the two pieces of T, 

f = fe(p) + fe(-p), (20) 

the first acting on the positive momentum subspace and the second acting on the negative momentum subspace, and 
combining them with a negative sign instead, 

f_ = T®(p) - f @(-p) . (21) 

This operator is self-adjoint, and has eigenstates |T) formed by the combination 

\T) = \T + ) + \-T-). (22) 

This time operator avoids Pauli's argument in a different manner. It is not conjugate to H but to an operator 
TL = sgn(p)H related to H by the change of the sign for negative momenta. Technically, because T is a maximal 
symmetric operator, TL cannot be its self-adjoint extension, and only gives the same result as T when acting on the 
subspace of positive momenta (there is a change of sign for the negative momentum subspace). 

So one can in principle construct self-adjoint variants of T. But what do they mean? and what new information 
can we extract from them? Unfortunately, neither the states \T ± ) nor the states T) transform according to Eq. 
(||). Equivalently, none of the self-adjoint operators discussed satisfies the covariance condition, and the "arrival 



time distributions" computed with them do not in general satisfy the basic physical requirement (12). To avoid this 
problem, in both approaches the domains of physical applicability of the self-adjoint operators have to be restricted 
with respect to the mathematical domains. For the approach by Grot, Rovelli and Tate the domain should be restricted 
to states with momentum support outside the regularization region around p — 0, (— e,e), so that in fact the arrival 
time distribution is again given by ([To|). Similarly, the approach associated with T_ is only physically meaningful 
for states with purely positive/negative momenta and the corresponding distributions are therefore contained in (|1C|). 
In summary, the apparent advantage of these two proposals is misleading, since in practice they are only applicable 
when their results are equivalent to the ones provided by the POVM related to the Aharonov-Bohm operator. In fact, 
the limitations on the domains of applicability of these two approaches are quite severe. The arrivals at a screen or 
detector will occur for all states in TL and not just for an especial set of states. A complete theory should provide the 
arrival time distribution in all cases. The important point is that the distribution associated with the POV measure 
applies for all states in TL (irrespectively of their behaviour at p — 0), so it is in this sense a more complete approach. 
This is perfectly compatible with D(T) ^ TL since the time operator is only one of its moments. Let us recall that a 
probability distribution exists independently of the existence of its moments. In fact in classical mechanics ensembles 
with non-zero probability at p — have no finite average arrival time, but the distribution is nevertheless well defined. 
In this respect only the POVM approach provides a correct classical limit. 

In summary, the distribution n(T) is satisfactory in many ways. There is however an important point that has 
not been considered yet. How do the eigenstates of T behave? Do they really represent states that arrive at a given 
time for a given position? The meaning of these eigenstates, although central, has not been sufficiently discussed. 
In particular, they have always been studied in momentum representation without paying attention to its coordinate 
representation and time dependence. These aspects are examined in the next section. Before doing so, it is noted that 
, because of the symmetry (|l6|), the coordinate representations of the time evolved states \T t ± } = exp(— iHt/h)\T ± } 
are related by 

(x\T t -) = (- x \T t + ). (23) 

so that it is enough to study one of the cases. Further simplification comes from the fact that, as a consequence of 
(|J), studying the time dependence of one of the states, say |T/ + ), from t = to the nominal arrival time t = T', 
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is equivalent to considering the sequence of eigenstates \T + ) from T = T' to T = 0. The space-time analysis that 
we shall carry out in the next section should be taken with some precaution because these are not square integrable 
states. As occurs with the continuum stationary states used in scattering theory, their physical interpretation requires 
the construction of normalizable wave packets peaked at one of them. 



III. COORDINATE REPRESENTATION OF THE EIGENSTATES 



The coordinate representation of the wavefunction (Q) is given by the integral, 

1/2 



<*l T + >=/ <*>Gfs) 



Jp 2 T/(2rnh) 



dp, 



(24) 



where delta function normalization, (x\p) = hr 1 ! 2 exp(ixp/h), is used for the plane waves. Let us assume for the time 
being that T > 0. The original path of integration is very inefficient numerically because of the rapid oscillations of 
the exponentials. The easiest way to calculate the integral is to deform the contour in the complex p-plane along the 
imaginary axis from the origin up to the intersection with the steepest descent path defined by 



pi = mx/T +p R , 



(25) 



(pr and pi are, respectively, the real and imaginary parts of p) and then to follow this steepest descent path rightwards 
to infinity in the first quadrant. The saddle is on the real axis at —mx/T. Two cases have to be distinguished: For 
x > the path does not cross the saddle and the value of the integral is small. In this case the origin is the only 
relevant critical point. Use of Watson's lemma pl| for large x and small T provides the leading term 



(x\T- 



2 .3/2 7r2 5/2 m l/2 ' 



OO . 



(26) 



note that this asymptotic behaviour is independent of T. But for x < the saddle becomes the dominant critical 
point. Retaining the leading term |l9[, 



(x\T + \ 



(m\x\/h)^ „ /4 



T 



(27) 



In the first case (x\T + ) decreases as a; -3 / 2 , whereas in the second one it increases as Ix] 1 ^ 2 asymptotically. In fact 
the integral can be expressed exactly in terms of parabolic cylinder functions, but the approximate critical point 
treatment just outlined is worthwhile since it allows a simple rationalization of the exact results. By means of the 
change of variable 



ap, 



where 



(p4|) becomes 



T 

mfi 



-1 1/2 
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1 
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where all square roots are calculated with a branch cut on the negative real axis, and 

/my 



xe 



-**/* (H!lY /2 



(28) 



(29) 



(30) 



(31) 



The integration path C in ( |30| ) goes now from to oo along the bisector of the second quadrant in the complex s 
plane. Note that z is the saddle point of exp(w) = exp(— s 2 /2 + sz). It lies on the bisector of the second/fourth 
quadrant for positive/negative x. The steepest descent paths from it are parallel to the real axis, and C lies on the 
border between "hill" and "valley" (wr = Real(w) = 0). It can however be deformed into a line just above the branch 
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cut since there is no contribution at oo. As the integrand across the branch cut simply changes sign, ( |30| ) is one half of 
the loop integral around the cut. In this fashion one of the integral forms of the parabolic cylinder function D^ 3 / 2 (z) 
can be recognized Q , 

^ T + ) = jJj^y /iD -^)- (32) 



The asymptotic behaviour for large \z\, see |21|], is in agreement with the expressions (|26|) and fl27|). The corresponding 
results for the case T < are simply obtained by using the symmetry 

(x\T + )* = (-x\-T + ). (33) 

The asymptotic behaviour for negative and positive x (growing and decaying respectively) changes abruptly at T = 0. 
Since \T+) does not represent a physical state vector, the discontinuity at T — is not problematic, but indicates again 
that a literal physical interpretation of these states is not allowed. We shall see later that normalized wavepackets 
formed with these states do not present this singularity. 

There is a region close to x — that cannot be described by the asymptotic formulae for large argument. In that 
region, however, the parabolic cylinder function can be expressed by means of a power series pOfl , 

r(-l/4) T(l/4) a 

D -V2(z) = nl/225/i - n l/2 2 3/4 Z + ( Z )■ ( 34 ) 

Combining ( |34| ) and (|3^ ) explicit expressions for (x = 0|T + ) and for its "flux" J(x — 0) can be obtained [p2fl . In 
particular, 

(27r) 3 T223/2 ■ l^J 

For T > 0, neither the wave function (x\T + ) nor the flux are zero at x = 0. Both the "density" \(x = 0\T + )| 2 and 
J(x = 0) grow monotonically as T — > 0. Figures la and lb illustrate all the dependences discussed. |(a;|T + }| 2 is 
depicted for a series of decreasing times T > as a function of x for two different scales: between x = —2 and x = 0.2 
(la), and between x = —0.2 and x = 0.2 (lb) (atomic units are used in all numerical examples). In the larger scale 
the density of the eigenstate is essentially a straight line, pivoting at x — 0, that approaches the vertical as T — > 0. 
The finer scale however shows that the arrival is not sharply defined. Even though, in a loose sense, "most of the 
wave" (x\T + ) passes from x < to x > at T = 0, there is a tail at x > present for an arbitrary T. Since these 
states are not normalized to one it is not possible to quantify the fraction of particles that can be found to the right 
of x = before T = 0. 

One might think that the eigenstates of the self-adjoint operators, which are orthogonal, could avoid this type of 
unsharpness |2j|. But they don't. For arbitrary values of e we have numerically checked (see the Appendix B) that 
there is a non vanishing density \(x\T + ) e | 2 at x > for fixed T. Similarly, the states (x\T) are also non zero on the 
right hand side. Using their defining equation (^2|), ( p3|) and (|33|), one finds 

(x\T) = 2Re{x\T+ > . (36) 

(It is to be noted that the corresponding time evolved state does not take this form. Instead, (x\T t ) = (x\e~ lHt ^ h \T) = 
(x\(T — t) + ) + (x\(T + 1) + )*.) Figure 2 represents the square of this quantity for a relatively large spatial interval 
(the small, but non-zero density for x > is not seen in this scale). The figure corresponding to the imaginary part 
is very similar. Note the increasingly rapid oscillation as x — > —oo. Cancellation of these oscillations occurs in linear 
combinations of the states (x\T) or (x\T + ) over a nonzero band leading to localization in a region closer to the 
origin, see Figure 3 and the discussion in the following section. There is a helpful classical association to understand 
the oscillatory pattern: high velocity particles (fast oscillations) have to start at longer distances and low velocity 
particles (slow oscillations) at shorter distances if they all have to arrive at X = simultaneously. 

IV. COORDINATE AND TIME DEPENDENCE OF THE NORMALIZED QUASI-EIGENSTATES 

We shall here construct normalized wavepackets by using a Gaussian distribution of the states (x\T[ + ) peaked at 
T' = TM. 



G 



/OO (T — T 1 ) 2 

e~ { x \Tl + )dT'. (37) 

-OO 

Carrying out the Gaussian integral over T' and then determining the normalization constant iV by evaluation of 
J Too d x I (^I^K^j T, AT)) | 2 gives the normalized wave packet 

o_l/4 1 A T\ 1/2 f°° 

(x\*(t;T,AT)) = — [ ;J / p i/2 e ^A e «(T- t )P 2 /2mR e -(AT)V/8m 2 r l 2 dp _ (38) 

This integral is readily evaluated numerically because of the exponential dependence on — p A . The momentum repre- 
sentation is given explicitly by 

<p|¥(t; T, AT)} = 27r 1 / 4 (Ar) 1 /2 ^ i/2 ^ (T _ t)p2/(2mri) ^ (Ar)V/(8m2?t2)e ^ 
(/im) 1 / 2 

As AT — > these states are orthogonal and satisfy the eigenvalue equation to any desired degree of accuracy: 

(mt;r,AT)Mt-,T,AT)) = w f^)={l (_^\ i T _ArjL ( 40 ) 



(p|T|*(0; T, AT)) = T(p|*(0; T, AT)) + C(AT) 5/2 , (41) 

where w(z) = exp(— z 2 )erfc( — iz) is the "w-function" pC|| . The moments (p 71 ) of the momentum distribution can be 
obtained from (|39|), 

«-(^)^(^)""' ! - 

Both the average energy (E) = {p 2 )/2m = h/(2n 3 / 2 AT) and the variance AT = [2" 1 - 7r" 1 ] 1 / 2 ^/AT diverge as AT 
goes to zero. 

Similarly, from (^8|) one finds 

(^-M^yv,,, (43) 

7r \mAl ) 

whereas higher moments diverge. (A finite "width" can however be defined as the half width at half height.) The 
average velocity {p/m}, which is also the velocity of the centroid (x) is given by r(3/4)[/i/(ATm7r 2 )] 1/ ' 2 . The behaviour 
of these states is arbitrarily close to what one could desire for an ideal arrival time eigenvector: They are normalized, 
and there is no discontinuity at t — T; for a fixed AT, they still obey a transformation law of the form (^J); the wave 
travels towards the origin with constant velocity, and the average position crosses the origin at time t = T, which is 
also the time when the spatial width of the wave packet attains its minimum value. Of course a certain unsharpness 
remains: From ( p^ ) one finds ( — x\^(2T — t;T, AT)) = (^(t;T, AT)\x), so the probability density has inversion 
symmetry in (x,t) with respect to the space time "origin" (0,T). In particular, for t = T, half the norm is to the 
right of x = independently of AT (The question of the existence of quantum states where the particle stays strictly 
on one side of X before T and on the other side after T is addressed in the Appendix C.) However, the wave packet 
density and flux are more and more peaked at the space-time point [x = 0, t = T) as AT — > 0, so that the passage 
of probability from left to right is sharp to any desired accuracy, 

* xv,9 f m \!/2r(3/8) 2 , % 

|<z = 0|*(i = T;T,AT)>| 2 =(— ) (44) 

r(3/8)r(5/8) 

J(x = 0,t = T)= Ar7r3/223/4 ■ (45) 

The normalized wavepackets constructed with the eigenstates ( p"7j ) have been examined in ]25| ], where it is reported 
that this arbitraryly sharp accuracy is not found in that case. 
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V. DISCUSSION 



Given the importance of the timed detection of particles at screens or in time-of- flight experiments, finding a theo- 
retical description of the arrival times seems imperative. Apparatus dependent results are available, for example via 
modelling the detection with complex absorbing potentials or other measurement models [p6|-p8[ , but it is reasonable 
to inquire if an intrinsic, apparatus-independent distribution can be naturally defined by means of the usual oper- 
ator approach to quantum mechanics. We have seen in sections II and III that the POVM and the corresponding 
distribution II (T) associated with the time operator T provide a rather satisfactory answer from the point of view 
of the properties satisfied: covariance with respect to time translations, minimum variance, appropriate symmetries, 
physically correct domain of applicability, and sharp space-time behaviour of the normalized quasi-eigenstates. 

One of the objections by Paul to T, which is also a shortcomming of the theories based on self-adjoint operators, 
namely the restrictive domain of the time operator(s) (which does not include, for example, states such as minimum 
uncertainty product Gaussians) is overcome by the POVM theory. An important point is to consider the primary 
object as the POVM (equivalently the resolution of the identity, or the arrival time distribution) rather than the 
operator. In this manner, the domain of applicability of the theory is TL, and the classical limit is correct. Moreover 
in this light the other problem indicated by Paul, namely the extreme sensibility of the expectation values of powers 
of T to small perturbations, is relatively unimportant. It is a fact that some quantities are very sensitive to certain 
small perturbations and we are simply dealing with one of them. The moments in the classical case would also suffer 
from such a sensitivity. The moral is that one should not pay as much attention to the moments of a time-of-arrival 
distribution (highly unstable with respect to small perturbations or changes in the apparatus resolution, and possibly 
divergent) but to general features of the distribution (peaks, global form, or width at half height for example). A 
consequence is that an uncertainty principle in terms of (T 2 ) is not of much use, since this quantity will generally 
diverge. If it does not, it will be too unstable with respect to small perturbations. It is preferable to express the 
uncertainty principle in terms of other measurements of width, such as the the half width at half height. 

Apart from the the positive features of n(T), it is also fair to point out several unclear points or open questions. For 
example, it is necessary to generalize the treatment to higher dimensions or to scattering problems. The connection 
between T1(T) and measurements (especially in non-classical cases) is also pending, see also |27|. We can however 
advance as a preliminary analysis two important features of a hypothetical measurement of n(T): The resolution of 
identity (||) implies the following structure of II(T), 

n(T) = |(T+|^ + )| 2 + |(T-|V,-)| 2 , (46) 

where \4>±) = Q(±p)\tp). This implies ignoring the "interference terms" \ip-)(ip+\ and so that in a hypo- 

thetical operational procedure to measure n(T), only the diagonal terms of the density operator contribute, 

Apart from the fact that the practical implementation of this diagonalization may be cumbersome, the neglect of 
interferences is not a desirable feature, since many different quantum states would give the same distribution. A 
second problematic aspect is associated with the interpretation of \(T + |V>+}| 2 as the contribution to II(T) from 
particles arriving from the left and of \{T — \ip-)\ 2 as the contribution from those arriving from the right. This is 
particularly evident when the backflow effect, 

J < for M = |V>+) , (48) 
J > for \ip) = IV--) , (49) 

occurs. Then II(T) assigns zero probability to arrivals from the "anomalous side", (e.g. from the right when \ip) = 
|V>+))- This implies that if \ip) = \ip + ) or \ip) — \tp — ) then either particles should be found to arrive either only from 
the left or only from the right, respectively, even during the time interval when J has the "wrong sign", or that the 
theory is appropriate when the "screen" is "one-sided" , failing to detect any particles arriving from the "anomalous 
side". Now consider the corresponding implication for the general state 1-0) = \ip + ) + \ip — ) with both \ip + ) and 
\ip — ) nonzero. One possibility is that the interference terms do not in fact contribute to the intrinsic arrival time 
distribution. The other is that the distribution ( (46] ) is only appropriate when the apparatus measures the sign of the 
momentum of each incident particle, thus collapsing the wavefunction of that particle either to \ip + ) or \ip-), and 
then switches on the appropriate one-sided detecting screen. 
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APPENDIX: A: LARGE T BEHAVIOUR OF U(T) 
Let us consider the asymptotic, large T behaviour of (T + \ip), 

1 /*oc 

(T+\<P) = J^Tj2j Q P 1/2 e-^ T ^ 2mh \p\iP)dp. (Al) 

Assumming that (p\ip) can be analytically continued into the fourth quadrant in the complex p-plane we deform the 
integral contour into the ray p — 7e~ 171 "/ 4 , (0 < 7 < 00). With the change g = r y 2 /(2mh) the integral takes the form 

{T + m = { -^^e-^ J™ dg e-^{g)9 VA , (A2) 

where the origin appears as the critical point. If, as g — > 0, ip(g) ~ cg u (where u is not necessarily an integer), use of 
Watson's lemma provides the dominant term, 

cr +W ~.-^ y r( £ff ) - (A3, 



APPENDIX: B: COORDINATE REPRESENTATION OF |T + } £ 

The integral over p for the coordinate representation {x\T + ) can be separated into two parts, from to s, I\, 
and from e to 00, J 2 . The second one can be obtained numerically with the treatment of section III, by deforming 
the contour and using the steepest descent path from the saddle. I\ can be expressed as, 



h = 



£ l-iTs 2 l(mK) pe 



f e ixp/h p (^~h) dp ( B i) 
Jo 



hm 1 / 2 

With the change of variable u = —ixp/% the incomplete gamma function is recognized, 

,3/2 



; ' = ^T 2 (W T(^> —ixe/h) (B2) 

£ 3/2 ~ {lxe/h) n 

hm i/2 Z_, n \(A + n) ' 



n=0 



where A = iTe 2 /(mh) + 1/2. It is easy to calculate an upper bound on the p < e contribution, I\{x,t;T), to 
< x\T t + >. The change of variable p = eu gives 

r, . ™ e3/2 f 1 du \.(txu e 2 tu 2 e 2 T, \1 mts 

hMT) = h^r 2 L ^ exp K— ~ 2^ + ^HJ- 

Hence, both the real and imaginary parts of h(x, t; T) are bounded in absolute value by 

hm 1 / 2 J u 1 / 2 hm 1 / 2 ( ' 

and \h(x,t;T)\ < {2ef/ 2 /{hm 1 / 2 ), independent of x,t and T. 
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APPENDIX: C: STRICTLY SHARP ARRIVAL STATES 



Are there quantum mechanical states where the particle stays strictly on one side of X before T and on the other 
side after T? We have seen, by examining the coordinate representation of the eigenstates of different time operators 
that none of the proposals (fulfilling cither the covariance condition or the self-adjointness condition) satisfies strictly 
this requirement. Is this an inherent limitation of standard quantum mechanics? We shall argue that there are no 
quantum states, pure or mixed, that satisfy fully this condition. To this end let us use the equivalent phase space 
Weyl-Wigner formalism. Its advantage for free motion dynamics is that the evolution kernel and dynamical equation 
of motion are equal in classical and quantum mechanics f?o|| . Therefore, a Liouville theorem applies, so that each 
phase space point carries its own "probabilistic weight" (that can be negative in the quantum case), so that in the 
intermediate calculations one may think and operate classically, the only difference being in the domains of states 
allowed in both mechanics and in the interpretation of the formalism ]30| , pT[ |. The Wigner distribution f(q,p) in the 
position-momentum phase space represents a valid quantum mechanical state if the associated density operator p is 
positive, see e.g. fl32| . Necessary conditions may be found for f(q,p) itself, such as 

\f(q,p)\ < (2/h) (CI) 

h [ f(q,p)dqdp<l. (C2) 



What kind of ensemble of classical particles, having negative positions at time t — 0, would arrive at the same time 
T at a point X = 0?. Since momentum is conserved, it is necessary a coordinated motion where faster particles start 
moving further away and slower particles start closer to X = so that they all arrive at the same time. The phase 
space density that satisfies these requisites is 

fT(x ,p;t = 0) = g(x )<d(-x )5{x +pT/m) , (C3) 



with g(xo) > 0. But this distribution is too singular to satisfy (CI) or (C2) 
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FIGURE CAPTIONS 

Figure la \(x\T + }| 2 versus x for T = 0.01 (long dashed line), 0.005 (short dashed line), and 0.001 (solid line); 
m = 1. All quantities in atomic units. 

Figure lb Same as Figure la for a smaller x interval. 
Figure 2 |(x|T}| 2 versus x for T = 0.01, m = 1. 

Figure 3 Probability density of normalized quasi-eigenstates, see (|37|), for T = 0.04, m = 1, AT = 0.002, t = 
(solid line), t — 0.02 (long dashed line), and t = 0.04 (short dashed line). 
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